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We relate the observable number of sources per solid angle and redshift to the underlying proper 
source density and velocity, background evolution and line-of-sight potentials. We give an exact 
result in the case of linearized perturbations assuming general relativity. This consistently includes 
contributions of the source density perturbations and redshift distortions, magnification, radial 
displacement, and various additional linear terms that are small on sub-horizon scales. In addition we 
calculate the effect on observed luminosities, and hence the result for sources observed as a function of 
flux, including magnification bias and radial-displacement effects. We give the corresponding linear 
■ result for a magnitude-limited survey at low redshift, and discuss the angular power spectrum of 

the total count distribution. We also calculate the cross-correlation with the CMB polarization and 
temperature including Doppler source terms, magnification, redshift distortions and other velocity 
effects for the sources, and discuss why the contribution of redshift distortions is generally small. 
Finally we relate the result for source number counts to that for the brightness of line radiation, for 
\& , example 21-cm radiation, from the sources. 



o 

u 



o 

H 



(N 
> 



I. INTRODUCTION 
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Redshift surveys of galaxies and other sources consist of a set of source positions along with measured redshifts and 
other information. With this information we can hope to learn much about structure formation and the background 
cosmology. However predictions for cosmological models are usually calculated as distributions over spatial hyper- 
surfaces. To relate to observation we have to map our light-cone observables into the underlying physical densities, 
velocities, and other relevant quantities. The aim of this paper is to derive exactly the relation between these quanti- 
ties, and provide a fully self-consistent calculation for scalar perturbations in the approximation of linearized general 
relativity. Relating the physical densities of sources to underlying cosmological parameters is then a separate further 
problem requiring detailed understanding of the source population and the bias. 
, This work arose from asking the question: when comparing observations of the matter power spectrum to theory, 
Q\ • in which gauge should the perturbations be calculated? For the case of current galaxy surveys, the answer is that 
fS| \ it does not really matter since, for observations well within our Hubble volume today, gauge effects are small (we 
observe wavenumbers with k/% 3> 1), and cosmic variance is large 1 . In practice, what is normally calculated is the 
| synchronous-gauge dark matter power spectrum since this is smooth in k as the Hubble scale is crossed, and then 
. some model is assumed for the bias between source counts and the dark matter density. However, it is well known 
that redshift distortions are important, which depend on the source velocities, so we are not really directly observing 
any single source power spectrum but a combination of various effects that contribute to the observed source counts 
' as a function of redshift. 

In this paper we calculate self-consistently (under simple assumptions about bias) all the various non-stochastic 
linear effects that enter calculations for the observed angular densities of sources per redshift. The most important of 
the additional effects have calculated separately before, but for completeness we also include all the velocity and post- 
Newtonian effects that appear in a full linear analysis. We also discuss the result for sources observed as a function 
of flux, and give a new result for a magnitude-limited survey including magnification bias, radial-displacement and 
source evolution effects. In addition we calculate the correlation with the cosmic microwave background temperature, 
including various velocity and post-Newtonian effects, as well as a new calculation of the cross-correlation with the 
CMB polarization. Doing a self-consistent analysis means we can be sure that no terms have been prematurely 
neglected (within the assumptions). In Appendix [a| we relate the result for source number counts to that for line 
radiation. 
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1 Note that the turnover in the matter power spectrum is determined roughly by the horizon size at matter-radiation equality; differences 
in gauge are important on scales approaching the horizon scale at the epoch of the observed sources, which is a much larger scale than 
the turnover for observed redshifts z <C 1000. 
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Differences between the general-relativistic and Newtonian predictions are in most cases beiow cosmic variance 
(depending on exactiy what is compared), though potentially a source of bias in cross-correlation studies. For the 
study of scale-dependent halo bias (as a probe of primordial non-Gaussianity jjj) using a joint analysis of multiple 
biased tracers, small corrections may be more important since there is only one underlying cosmological perturbation, 
and hence in principle no cosmic variance on the difference between source counts |^]-^| . In particular there is a danger 
that an incorrect calculation of the large-scale power spectrum could give spurious evidence for scale-dependent bias 
due to the neglect of various non-negligible effects that are present in the full analysis. 

Previous work has examined perturbations to the luminosity distance in detail Reference pif considers the 

total source count to up to a certain redshift in a specific cosmological model, but does not give results for counts 
as a function of flux or magnitude limit. More recently, while this work was in progress, Refs ]TT|-|l3[ have given an 
analysis more similar to this paper, though our presentation is rather different and our results do not all agree. In 
particular our numerical results are quantitatively different, in part due to our more consistent treatment of bias. Our 
numerical "CAMB sources" code for calculating source count angular power spectra and cross-correlations with the 
CMB and galaxy weak lensing is publicly available 2 . 

The outline of this paper is as follows. In SeoO we give a general non-perturbative result for the observed source 
densities as a function of redshift, then in Sec. HI we calculate the result for linear scalar perturbations. We use the 
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Newtonian gauge, however since we are calculating an observable the final result is gauge invariant. In Sec. IV 
include the source luminosities and derive the observed source density for a magnitude-limited survey. In Sec.^ we 
discuss the dominant terms in the full result, explain how we perform numerical calculations over a source window 
function, and present some typical numerical results. We also discuss evolving source populations, and the distinction 



between the selection function and the underlying background source redshift distribution. Finally in Sec. VI we give 
results for the cross-correlation of source counts with the CMB temperature and polarization, including all the various 
effects that can be important for both low and high redshift sources. In Appendix ^ we relate the result for source 
number counts to that for line radiation (e.g. 21-cm radiation). 

II. PROJECTED NUMBER COUNTS 

We first calculate the observable angular source density count as a function of redshift and angle on the sky in terms 
of the physical quantities governing the source density, velocity and line-of-sight evolution. Readers not interested in 



details of the derivation can skip to the main linearized result in Eq. (30) or ( p7| ) below. Throughout we use units 
with c = 1. Our derivation below, up to and including Eq. (p^), holds for any metric theory of gravity. However, all 
numerical results assume general relativity. 

We consider a general population of objects with current 4- vector J a = n s u® where n s is the proper source number 
density (i.e. in the rest space of the source 4- velocity m°). These are observed in projection by an observer at A who 
has 4-velocity u a oA . The observer records n(n, z)dzdfl A, the number of sources in direction n over solid angle dCl A 
with redshift z in a range dz. Let k a = dx a /dA be the wave 4- vector of a light ray on the past lightcone through A. If 
the ray intersects a source at affine parameter A, the redshift of the source is determined from the observed frequency 
and the emitted frequency in the source rest-frame by 

1 + z= (*«<)U (1) 

(kaU^ A )\\ A 

The infinitesimal angular separation 80 j of two rays in a bundle at A, and the ray-orthogonal connecting vector 
between them at affine parameter A along the rays, are related by the Jacobi map T>u(X) (see Refs Hj-Q): 

£i(A) = T>ij(X)56j, (2) 

where the indices here are components in a 2-dimensional orthonormal basis {Ef} orthogonal to k a . The Jacobi map 
is determined by its evolution equation along the ray, 



TikVkj, (3) 



dA 2 

where the optical tidal matrix is defined in terms of the Riemann tensor by 7},/ = — £7j J E7^fc CL fc d J? cl b cc ;. 



2 lttp : //camb. info/sources/ 
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We can now use the Jacobi map to relate observed ray angles to (invariant) areas on the wavefront: a ray bundle 
with solid angle dtt A (as seen by u^ A ) has invariant area det£> df2 0j 4, where det£> D is the determinant of the Jacobi 
map in the u" A frame. In incrementing by dA, the wavefront advances by a proper distance d\k a u® in the source 
rest-frame, and so sweeps up n s detT> k a u^dXdn o A sources. The general result then follows: 



n(n,z) = deW k a J a 



dA 

dz 



(4) 



where the right-hand side is evaluated at redshift z along the line of sight with direction n at the observer. 

Under a change of observer, u a oA — > u a oA (so that n — > n and z — > z), conservation of the number of sources gives 
n(n, z)d£l Adz = n(n, z)dil a Adz. It follows that 

n(£, z) = (^A n(n, z) = 7 3 (1 + n • v rcl ) 3 n(n, z), (5) 

V KaU oA J 

where v rc i is the relative velocity of the observers and 7 is the associated Lorentz factor. 

The result for number counts is related in detail to that for the brightness of diffuse line radiation from the same 
sources in Appendix [A|; in summary, the fraction of photons received from a source scales with 1 /det£> (inverse-square 
law in the background), so the sky brightness is simply oc n(n, z)/det£> . 



III. PROJECTED NUMBER COUNTS IN FLAT, ALMOST-FRW MODELS 

We aim to evaluate Eq. (|j) to linear order for a fiat, almost-Friedmann-Robertson- Walker (FRW) model with scalar 
perturbations. For convenience, we work in the conformal-Newtonian gauge with metric 

ds 2 = a 2 { v )[(l + 2V>)<V - (1 - 2<f>)5 ij dx i dx j }. (6) 

We take a zero-shear velocity field u a along d v so that w M = a _1 (l — iP)6q and it M = a(l + - 0)<5 A ,o- This velocity field 
is the zeroth clement of an orthonormal tetrad which we take to be (X ) a = u a and X { = cr^l + <f)di. 

Decomposing the wavevector k a = dx a /dX into a direction e a and frequency k ■ u = e/a relative to u a , we have 

^ = (l + + ^)e, ^L=a- 2 e(l-^), (7) 

where the three- vector e comprises the spatial components of the propagation direction on the spatial triad Xi. The 
geodesic equation reduces to a simple propagation equation for e, 

de 

^ = -V ± (<P + tP), (8) 
dr/ 

where = V — ee ■ V, and an equation for the evolution of the comoving frequency: 

-=-e/ + e(^ + # (9) 
dr; dr] 

Combining Eqs (0) and (||) and solving for the photon path gives 

rn rv 

x(n;r ? ) = -e A (r ?A -r/)+e A / (cf> + <ip) dr}' - I (r? - »/)V±(0 + V)<V- (10) 

Jr/A Ji] A 

Note that at the observation point dx/dr/ = (1 + <f>+ ip) A eA, so that — eA is the line-of-sight direction for a Newtonian- 
gauge observer (u A ). The second term on the right of Eq. ( |l0| ) is a radial displacement and corresponds to the usual 
(Shapiro) time delay. The third term is the usual transverse (lensing) displacement. 

We shall need the radial displacement and perturbation to the conformal time at redshift z along the line of sight. 
Writing the source 4-velocity as it" = u a + v a , we have u£ = a _1 [l — "07^] where v % are the orthonormal-triad 
components of v a , and similarly for the observer's 4-velocity = u a + d°. It follows that the observed redshift of a 
source is 

i + z = ^i_L(i + fi.[ v _ VoA ]) i (11) 

a eA 
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where n = — e a A is the line-of-sight direction on the observer's triad formed from Lorentz boosting the (Xo) a and the 
(Xi) a . At zeroth order, b a = e^. The ratio of energies follows from integrating Eq. (^): 

— = 1 + i>A - i> + [\<P + ip) drf, (12) 



(A 



>]A 



which has the usual Sachs- Wolfe and integrated Sachs- Wolfe (ISW) contributions. The redshift at rj along the line of 
sight is therefore 

1 + z = ( 1 + i> A - i> + f + 1>) dr]' + n • [v - v oA ] J . (13) 



a(rj) 



Setting rj — rj* + 5rj for a source at observed redshift z*, where 1 + z* = 0^/0(77*), we must have 

f 7 >* . 

H(r}*)6r} = ip A - 1/> + (0 + V>)dr/ + n- [v - v„a], (14) 

Jr)A 

where H is the conformal Hubble parameter and the terms on the right are evaluated on the zero-order lightcone at 
position + h.(r)A — f)*) at time 77*. The radial position of a photon at observed source redshift z* then follows from 
Eq.©: 

X(n,z*) = X* +$X 

= VA - V* - - / {4> + ip) drf . (15) 

Jr\A 

To determine the perturbed value of det£> = det'D&Q a / d£l A^ we note that the Jacobi map is symmetric for linear 
scalar perturbations. It thus takes the form 

»»- ( D/ !; 71 ^ 71 ). ™ 

where V is the trace of the map and 71 and 72 are the (first-order) components of the shear. Evaluating the 
determinant, we find at linear order that detl? = (V/2) 2 . The trace is given in terms of the convergence K by (see 
e.g. Ref. g§) 

V{h, V )/2 = X (n, V )a( V )[l - <j> - K {n,r])], (17) 
where x is the perturbed radial position at time 77 and 3 

2«(n,»,) = -V| f ^T^ A ^ + ^drj'. (18) 

The Laplacian here is on the unit sphere. We actually require detl? at given observed redshift] this is given by 

det£>(n,z«) =a 2 ( V *) X l U + 2^ + 2H5r) - 20 - 2k\ , (19) 

and hence 



det2? (n, z*) = dct2?(n, 2*) 



dfl a/ 



= a 2 (v*)xl ( X + 2 ~ + 2H67 1 ~ 20 - 2k + 2n • v oA ^j . (20) 

We also require the perturbed value of k a J a — n s k a u^ = n s (l + z)e A (l + n • v a)/o,a where, recall, n s is the proper 
(physical rather than comoving) number density of sources in their rest frame. Expanding n s into a background part, 



Note that other authors (e.g. [[l2| [l?]]) have defined the convergence differently so that it includes additional terms that we are including 
separately. 
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n 8 , and a perturbed part n s 8 n , and evaluating the background part at the perturbed time for a source at redshift z*, 
we have 



fc a J Q (n, z*) 



(1 + z*)eAn s (rj* 



1 H — -£77 + S n + n • v, 



oA 



(21) 



The final term in Eq. (Q) is dA/dz. We can evaluate this by differentiating Eq. (filf): noting that 

a 2 

■(l-ip-n-[v- v a]), 



d??/dA = e(l-7/0/a 2 

6a(1 + «) 



aaA 



we find 



dA 
dl 



(n,z„) 



a(r]^)aA 



e A H(r]*)(l + z*f 



U 

n 



H\6ri 



1 di/> 1 



1 , dv 



H dr; H 



+ —(4> + "0) + ttA ■ ~\ h "0 + n ■ [v - v 



H dry 



Putting these results together, the background result for the number counts per solid angle and redshift is 4 



n(z) = (ax) 



2 a 2 n s 



X 3 - 

?a n s , 



a A n H{l + z) 

where all quantities are evaluated at the source redshift z, with fractional perturbation 



(22) 



.(23) 



(24) 



Ti ■ oy 

A„(n, z) = 5 n + —Si] + 2— + 2H5r] - 2(f) - 2k + 3n ■ v oA 
n s X 



H 
H 



H\6ri 



1 dip 1 

h cb7 + n 



(0 + 



1 „ dv 
H n 'd^ 



■ V> + n ■ [v - v oA ]- 



(25) 



For conserved sources n s /n s = — 3%, so the terms due to constant expansion cancel. More generally, we can write 
A„(n, z) in terms of the background comoving density (a 3 n s ) evolution 

a /» r 1 » dv £y „ dln(a 3 n s ) . 
A n n, z) = <J„ + -n • — + 2^ - 2 K + y - s -±5r, 

+ 3n • v oA + n • [v - v oA ] - —Srj - + ^(0 + ip) + <ip - 20, (26) 

where terms have been ordered roughly in order of importance on sub-Hubble scales. This result holds for any metric 
theory of gravity. However, if we assume the general-relativistic velocity evolution equation for non-interacting cold 
particles, 



v + Kv + V?/> = 0, 



we can simplify further to give 



1 dv 5x 

A„(n, z) = 5 n - — n • — + 2 2k - 

H dx X 



dln(a 3 n s ) H 



Srj + 2n- v oA + —<f> + ip - 2(/>. 
ft 



(27) 



(28) 



The considerable simplification here is from a cancellation between the change in the Doppler shift across the source 
volume element due to the time evolution of n ■ v and the change in the Sachs- Wolfe effect due to the radial gradient 
in ip. 



4 This is simply the product of the comoving source number density (a 3 n s ) and the differential comoving volume element per redshift 
and solid angle. 
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FIG. 1: Three-dimensional power spectrum of the Newtonian-gauge counts source 5 n (solid) and the synchronous gauge source 
b5^ n (dashed) at z = and z — 3 with bias b — 1.5 (relative to the matter perturbation 5„ n ) assuming constant comoving 
source density. The dash-dotted line shows the spectrum for S n re-scaled (incorrectly) from z — 3 to z — with the scale- 
independent growth factor D(rj) appropriate to <5„ n . On large scales, S n grows more slowly in time than S^ n due to the velocity 
term in Eq. (Bit) which goes as HD. Note that neither power spectrum is directly observable, but the difference between them 
illustrates the importance of using a full analysis on scales k <C O.OlMpc -1 . 



There is a dipole contribution to A„(n, z) from the observer's velocity v a- Isolating these terms (including those 
implicit in Srj and 5x), their contribution can be written as 



(n, z) = 3n • w oA - n • w oA 



1 dlnn 
H dz 



(29) 



This is consistent with the linearized form of Eq. (g): the first term on the right is from (1 + n ■ v rc i.) 3 while the 
second term is from correcting the observed redshift to undo the effect of the observer's velocity when evaluating the 
background source counts n(z). Assuming the non-kinematic CMB dipole is O(10 -5 ), the local peculiar velocity w oA 
can be determined to 0(1O~ 5 ) from measurements of the observed (kinematic plus primordial) CMB dipole. If we 
boost to this frame the n • w q a terms then vanish. 

Expanding Eq. (Eq) and dropping local (monopole and dipole) terms the result is 



* /a i r 1 » dv 

A„(n, z) = 5 n - — n • 2k 

H dx 



2 dln(a 3 n s ) H 



-1>- 



HA 



((f) + ip)dr] - 
2 fVA 
X 



n • v 



^)&r,+ ^ + i)-2(l>. (30) 



The terms in this equation have a simple interpretation. The first is the perturbation to the source number density. 
The second is the usual redshift-space distortion due to the effect of source velocities on the differential volume 
element. The convergence term arises from the transverse differential effect of gravitational lensing on the volume 
element while the following two terms are from evaluating the zero-order result, Eq. (^4|), at the perturbed position 
and time appropriate to the observed redshift (including the radial shift due to lensing). The remaining potential 
terms are from the perturbation to the spatial metric {ip) in converting dz to a proper time interval at source, the 
differential of that part of the ISW effect across the source volume {(j>) that is not cancelled by velocity evolution (see 
above), and the effect of the spatial metric at source {2(f) on transverse distances there. 

The source number perturbation S n is the Newtonian-gauge quantity. In the simplest case we might expect the 
synchronous-gauge number perturbation S 1 ^ 11 to be related to the matter density perturbation via a constant bias 
b, <5^ yn = b5^ n , since it is the comoving perturbation (synchronous gauge when perturbations are dominated by 
pressure-free matter) that enters directly into the Poisson equation Jl8, 19 1. In this case the Newtonian-gauge number 
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density perturbation is, in Fourier space 5 , 



5n = b6%* + 



d ln(a 3 n s 
dij 



v 

k' 



(31) 



where v is the Newtonian-gauge velocity of the sources, assumed to follow the matter velocity (no velocity bias). The 
difference between the Newtonian- and synchronous- gauge matter power spectra is shown in Fig. [I], and illustrates 
that it is important to make the distinction on scales approaching the horizon scale. 

Our prescription for the bias differs from that in the papers by Yoo et al. [[12], where the assumption is made 
that the number density of sources n s — F(p m ) is some function of the local matter density. At the spacetimc 
point corresponding to the perturbed position of a source observed at redshift z, they linearise this relation about 
the background matter density at redshift z to define a linear bias relation. This singles out the zero-redshift- 
perturbation gauge and gives the local relation S n = bS m in that gauge. However, it seems clear that the bias is 
physically due to the local physics of structure formation, and nothing to do with how we observe it: hence in as much 
that bias is a reasonable model at all, its definition cannot depend on the zero-redshift-perturbation observational 
gauge. Note further that the synchronous gauge density is also what is calculated in Newtonian collapse and TV-body 
simulations |l|, |2l], ^| . For a recent discussion of bias in the general-relativistic context, see p3f . 

The fractional perturbation in the counts is not directly observable since we cannot observe n(z). Rather, we 
measure anisotropy in the counts via n(n, z)/(n(n, z)}„ — 1 = A„(n, z) — (A„(n, where denotes the angular 
average over the survey area. With full-sky coverage, the I > multipoles of the observed anisotropy equal those of 
A„(n, z). More generally, the difference in the observed fractional counts between two directions equals the difference 
in A n (n, z). Note also that the I > 1 moments of A„(n, z) are independent of the observer's velocity in linear theory. 



IV. LUMINOSITIES 



The linear-theory results derived in the previous section can also be applied to calculate observables other than 
total source number counts. Of particular interest here is the luminosity distance: how it affects the flux received 
from each source, and the predictions for the observed angular densities of a magnitude-limited source sample. 

If a source with 4-velocity a" is known to radiate an energy dE = L s dr s isotropically in its rest frame in 
proper time dr s in the form of photons with energy (k a Ug)\\, the number emitted within solid angle dfi s is 
dN = L s dT s dSl s /(4:Tr(k a Us)\\). These are collected by an observer at A, with 4-velocity u^ A and equipped with 
a detector of collecting area dA A , in time dr Q A with energy {k a ut & ) I \ ■ Products of areas and solid angles are related 
by reciprocity, (1 + z) 2 dQ s dA s — dfl AdAA (see e.g. Ref. pBj), and the Jacobi map relates areas and angles via 
dA s = detT> dQ A- The energy received at A is then given by 



dE oA = (k a u a )\ XA dN 



L s dT oA dA A 
4tt(1 + z) 4 det2V 



(32) 



where we have used dT A/dr s = 1 + z. 
the luminosity distance, given by 



If there are well-understood sources with known L s we can therefore measure 



d L (h,z) = (l + z) 2 {detV } 1/2 , 



(33) 



for each source. Our results for the linear-theory Jacobi map at given observed redshift could then be used straight- 
forwardly to work out the observed luminosity distance in a perturbed universe; see Refs 

The source number counts can also be measured as a function of source flux, F = dE oA j (dA A dT oA ) , so an observer 
u oA sees n(h, z,lnF)dil Adz sources per logarithmic flux interval dlnF. If the current of sources per log luminosity 
L s is J a (\nL s ) = Ugn s (x,r]AnL s ), we have the observed density 



n(n,z,lnF) = deW k a J a (In L s ) 



dX 

dz 



where 



L, = 4ttF(1 



z) 4 detV D . 



(34) 



(35) 



Our Fourier conventions follow Ref. |20 
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For populations that have smooth distributions in source luminosity, the linear perturbed result is then the same 
as before (Eq. ^) , with some extra terms depending on the slope of the source luminosity function: 



A„(n,z,lnF) = A„(n, z)\ laLs 
= A„(n,z)| ln£s 



dhin s 



d\i\L s 
31nn s 



91nL s 



lnL 3 



lnL 3 



detD 

2 s_x 

X.* 



1 



2H5ri - 2(f) - 2k + 2n • v oA 



(36) 



where A„(n, z)| ln £ a is evaluated for sources with luminosity L s = 4irF(l + z) a 2 (?y*)x*j i- e - the luminosity corre- 
sponding to flux F in the background. 

Apparent magnitudes are defined with conventional factors so that m = — 2.51og 10 i 7 ' + const. For a magnitude- 
limited survey observing all N(ii,z,m < m*)dzdfl A sources with magnitudes m < to* over redshift interval dz and 
in solid angle dil A, the full result for the fractional perturbation to N, neglecting only local terms, is 



1 dv 

A N (h, z, to < m«) = S N (L > L s *) - — n • — + (5s - 2) 

n dx 



1 

K 

X 



VA 



+ 



2 -5s dln[a 3 N{L > L s ,)] U_ 



>!A 



((f) + tp)drj 
ip)drj — n ■ v 



+ ^ + V + (5s-2)^. (37) 



Here, N(j],L > L s *) is the background number density of sources with luminosity exceeding L s * and its fractional 
perturbation is denoted i5at(x, r],L > L s *). The quantity 



s(z, m*) 



_ d log 10 N(z, m < m„) 
dm* 



n s (r],\aL s ») 
2.hN(mL > L„Y 



(38) 



where the terms in the final expression are evaluated at the background value of ?y corresponding to z. Note that we 
have assumed no line-of-sight scattering. 



NUMERICAL AND APPROXIMATE RESULTS 



For non-evolving source populations, Eq. (37) can be approximated on small scales as 



Ajv(n, z,m < m») 



)N 



1 . 9v 

n n 'dx~ 



n- [v - Vqa\ 



(2 - 5s). 



(39) 



The first term is the obvious underlying perturbation and the second the usual redshift distortion due to velocity 
gradients [ ^|, ^5) . The third term has two contributions from lensing: the first from the effect of lensing magnification 
of the background sources ]26|], and the second slope-dependent correction is the well-known magnification bias f27| - 
|30| . The first of the two contributions proportional to l/("Hx) comes from the change in the number of background 
sources per observed solid angle as radial distances are changed [^4|, |3l|-|3^|. The second slope-dependent term is 
the effect of source-dimming due to apparent radial displacements caused by source velocities ||: a source moving 
towards us is further away than it appears from its redshift and hence has to be more luminous to give the same 
observed flux. These 0(v/Hx) terms may be non-negligible compared to the neglected terms if the sources are close, 
X < KT 1 (roughly z < 1). 

For numerical work results can be Fourier transformed and expanded into multipoles in k and I space. For example, 
Eq. (p9]) gives 



A NJ (k,z) « 5 N ji{k X ) + ^-ji"(k x ) + (2 - 5s) 



fi(kx) 



XX 



(40) 



for I > 1. The full result can be calculated similarly, but for brevity we do not quote it here. Integrating over an 
observed redshift window function, W(z), gives a total fractional perturbation in that window 



A£,(fc)=/ dzW(z)A Nil (k,z)= dTjWfaJAjv,,^,*,), 
Jo Jo 



>IA 



(41) 
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FIG. 2: Angular power spectra of total number counts, based on the approximation of Eq. (|39|), for Gaussian window functions 
W(z) at z = 0.45 (left; a z = 0.03 and constant bias b — 1.5) and z = 0.6 (centre; o z = 0.05 and b — 1.95), and their 
cross-correlation (right) using non-linear corrections from Halofit [Q. Solid lines have no lensing (2 — 5s = 0), dashed lines 
have lensing but no magnification bias (s = 0), and dot-dashed lines have s = 0.6. Note that these are barely distinguishable 
in the auto-power spectra. The window functions are similar to those actually measured using LRG surveys |37[ |38| . Lensing 
is a significant source of correlation when the cross-correlation is otherwise small, c.f. Ref. pq ]. The effect is well below cosmic 
variance on an individual I, but above cosmic variance over a range A; ~ 100. For s — the effect of magnification is negative 
since a magnified area has less sources per solid angle. 



where z„ is the redshift at conformal time r\ in the background and W(rf) = (1 + z)HW(z). Note that here W(z) is 
the normalized total distribution of counts in the window, rather than a probability distribution for observing sources; 
if the number of sources and selection function were constant then W(z) oc fi{z) = a 4 x 2 n s /H| z oc ax 2 /"H\z- Defining 



2 -5s dln[a 3 N(L > L s 
1 5s 



no?] 



2L 

u 2 



(42) 



J n 



and switching some integration orders, the (full) window-integrated counts transfer function can be calculated using 

r V a 



^nM = I d V 
Jo 



in '/) I s Njl (k x ) + ^ii"(*x)) + w Sv (v) WMkx) + + + ^M^x) f Q w 6v W)drf 



, )+ , )M/ , x) ,/>- 5s )^ + ^ 



x__?c, 2 _ 5 s W ,,, drj , 

o XX 



1 • 



+ W( V ) Jl (k x ) ( -0 + ^ + (5s -2)0 



. (43) 



where \' = tja — f]' (and similarly for \) . We can further integrate by parts to obtain an integral of a source against 
jl(kx) (assuming W(z) goes to zero at both ends). The linear-theory angular power spectrum can then be calculated 
using standard line-of-sight Boltzmann codes, giving angular power spectra 



Cf w ' =4tt / d]nW> w (fc)A#,(fc)A!5r,(fc) 



, w' , 



(44) 



for the cross-correlation between counts in windows W{z) and W'(z). Here Vfi(k) is the dimensionless power spectrum 
of the primordial curvature perturbation and the transfer function A^ ; (fc) is for unit initial curvature perturbation. 
Cross correlations with the CMB, weak lensing, 21-cm or other sources can be calculated similarly. Corrections for 
non-linear evolution can be accounted for approximately by an appropriate re-scaling inside the k integral |35| . 

We first consider the effects of the magnification terms in the approximate result, Eq. (|3^). These are generally small, 
but do increase the correlation of counts on different redshift slices. Figure ^ shows that the effect can be significant 
for widely separated slices on small scales when the correlation is otherwise small. As discussed by Ref. JM the 
magnification also has a small effect on the cross-correlation with the CMB (see later discussion in Sec. |VI[ ). For 
sources with bias b the relative importance of the magnification terms depends on (2 — 5s) jb. Examples for various 
surveys are given in Ref. p9[ , and an example of the large-scale effect is shown in Fig. |^ when 2 — 5s = 1. 
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FIG. 3: Angular power spectra for all sources in a Gaussian window function at z — 0.1 (a z — 0.01, b — 1), with (thick 
solid) and without (dashed) the 0(n • v/Hx) radial-displacement term. The dot-dashed line shows the equivalent result for 
a magnitude- limited survey with constant 2 — 5s = 1, and the bottom panel shows the fractional differences compared to the 
result with no radial-displacement terms. 



Now consider the effect of the 25x1 X ~ — 2n • [v — v a]/ (Hx) radial-displacement term. This is generally small for 
distant sources, but can be more important when the sources at are low redshift, as discussed in Refs J3l|— 133|] . The 
effect on the large scale power spectrum is shown in Fig. |3[ Having a magnitude limit with small positive s partly 
cancels the low-redshift radial-displacement effect: a radial displacement gives a larger volume per solid angle, but this 
is compensated by the fact that the sources appear dimmer because they are further away, so less sources are seen. The 
effect on the power spectrum for single redshift slices is well below cosmic variance, but may be important for a full 
redshift survey. The effect of radial displacements on the all-sources correlation function is discussed in Refs [j32|, 33| . 
If there are other selection effects in addition to the magnitude, for example on galaxy size |40 , orientation El or 
dust extinction along the line of sight [Q, they should also be included. Compared to the radial displacement, the 
effect of the other n • [v — v q a] terms on the power spectrum is down by a factor of 0(Hx) 2 ~ 60 at redshift z = 0.15 
(assuming the source population is not rapidly evolving). At z > 1, the n ■ [v — v a] terms are of equivalent size and 
all small. 

Figure [| shows the fractional errors that can be made if various terms are neglected when calculating angular power 
spectra for relatively broad redshift window functions. General-relativistic potential terms are only significant at very 
low multipoles for high redshift (i.e. where the contribution of large-scale modes approaching the Hubble scale can be 
important). Radial displacement terms are important at low redshift, and lensing (being cumulative) at moderate and 
high redshift. Although the post-Newtonian effects are generally small compared to cosmic variance, the additional 
terms entering the observed counts of biased density tracers is independent of bias. If multiple populations tracing 
the same density field with different bias are being used in a joint analysis the cosmic variance only enters into one 
eigenvalue of the sample covariance and hence small theory errors can potentially be rather more important. 



A. Total counts distribution 



We can calculate the angular power spectrum for any combination of redshift window functions W(z). In general 
these windows are not directly related to the underlying physical source distribution n(z), either because we have 
chosen to divide up the data into different redshift bins, or because of observational (selection) issues. However we 
can also consider the special case where we observe all of the sources, and hence wish to calculate the total angular 
counts power spectrum over the entire distribution, W(z) oc n{z) — a 4 x 2 n s /TL\ z - 
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FIG. 4: Fractional error compared to the full result for the counts angular power spectrum Ci with broad Gaussian window 
functions peaking at various redshifts with cr z = 0.32. The error shown is that obtained when various types of term in the full 
result ofEq. @ are individually neglected. The 'redshift' term is for redshift-distortions (radial derivative of velocity), 'velocity' 
terms are proportional to n ■ v and include the radial displacement effect, 'lensing' is the convergence term, and 'potentials' 
includes the effects of gravitational potentials at the source, time delay and the ISW. Solid and dashed lines indicate terms whose 
neglect reduces and increases the power spectrum respectively. The 'standard' result is the approximation given by Eq. dis| ) 
and for this the fractional error compared to the full result is shown (with solid/dashed lines denoting an excess/deficit). There 
is no source evolution and (unrealistically) b — 1 and s = 0. 



Assuming s — 0, the velocity terms in Eq. (|30| ) can be written 

1 „ d\ d\n[x 2 a 3 h s /n} n • v 

A « (n ' z) = -n n -T x + W 

1 „ <9v dln[(l + z)fi(z)] n ■ v 

= n • — — I ' x n . (45) 

n d X dr, n v ' 

Noting that dv/di] = v — dw jd\ and dz = — (1 + z)%drj in the background, the perturbation to the velocity terms in 
the total counts become 

dzn(z)A v n (h,z) = drj (J^ii 1 + z)n(z)n ■ v] - (1 + z)n(z)h ■ vj . (46) 
Integrating the total derivative and noting that the boundary terms vanish, we find 

dzn(z)Al(h,z) = - / dr; (1 + z)n(z)h ■ v. (47) 
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FIG. 5: The fractional difference in the counts angular power spectrum (left) and CMB temperature cross-correlation (right), 
compared with the 'standard' result calculated using Eq. ( pl9| ) for the distribution given in Eq. (^8|) . Thick lines are our limiting 
full results for the case of no source evolution (solid), in which case the source selection function is non-trivial, and the case 
when all sources are observed (dashed), so that P(z) = n(z) implying source evolution. The thin red line shows the change 
in the standard result from dropping the redshift-distortion term; this modification to th e 'standard' result brings it close to 
the full result when all sources are observed as expected from the arguments in Sec. V A. The blue dashed-dot line shows the 
change in the standard result when S s ^ n is replaced with the Newtonian gauge Sn (no source evolution), which is significantly 
larger than the difference between the standard and full results. 



Since all the sources are being observed, the angular number density is not affected by changes in the apparent radial 
distance of each source; the only remaining contributions come from time evolution between the light cone and the 
perturbed positions. The result is generally small: a fractional perturbation of 0(v) rather than <D(kv/H) from the 
redshift-distortion term alone. This also means that when calculating the distribution for all sources it is important 
to include not only the dv/dx term, but all the source velocity terms since these nearly cancel. Alternatively in 
many cases the velocity terms can be neglected entirely, though near cancellation of the velocity terms is a useful 
numerical consistency check. Similar comments apply to the potential contributions to the change in redshift distance. 
However, on the largest scales it is still important to use the correct Newtonian gauge density source, 8 n , rather than 
approximating it as proportional to the synchronous-gauge density perturbation. 



B. Selection function or count distribution? 



Given some observed angle- averaged count distribution P(z), the predicted angular power spectrum for the sample 
depends on the interpretation of P(z): if P(z) — n(z) (i.e. we are observing all sources), then we are calculating 
the total counts distribution; on the other hand the shape of P(z) may be mainly due to the observational selection 
function, in which case the underlying total physical count distribution n{z) then needs to be specified separately 
if the source-evolution terms in Eqs (^0|) and ( |37| ) are to be included correctly. In general there will be both an 
observational selection function and some change to the source populations with redshift, and hence P(z) alone does 
not fully specify the problem. 

As an example, we follow Ref. |12]] by considering a source distribution for a photometric quasar sample, 



P(z) oc z a exp 



(48) 



with (a, (3, zq) = (3, 13, 3.358) and hence mean and the peak redshifts of 2.66 and 3.0 respectively. We take constant 
b = 2 and 2 — 5s = —0.1. We compare the full numerical results for the extreme cases of P(z) = n(z) oc W(z) (i.e. 
observe all sources) and the limit in which the source population is not evolving (a 3 fi s = const., so P(z) — n(z)p(z) 
where p(z) is the probability of including a given source at redshift z and, from Eq. (53), n(z) oc ax 2 /H\z)- For 
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reference we compare to the 'standard' result (with no source evolution) 



A« d s 6%" - -a ■ - (1 - 5s)k = bS%° - -n ■ - (2 - 5»)«. (49) 



Our results, presented in Fig. |J can be compared with Fig. 2 of Ref. |12| , and show a significantly smaller difference 
between the full result and the 'standard' result for the equivalent cases. The main difference is our alternative 
treatment of bias (Eq. |3l]) and the explicit treatment (or not) of source evolution terms (c.f. Refs JTTJ, ^2j). Nonetheless, 
we find that the difference between the full and standard result is still several percent on the largest scales, and larger 
for the CMB cross-correlation where there are additional contributions from correlations between count velocity terms 



and the CMB reionization Doppler signal, and large-mode early-time Sachs- Wolfe contributions (see Sec. VI). It is 
also clearly important that the redshift distortion/source evolution terms are incorporated in a consistent manner or 
much larger differences can be obtained. 

VI. CMB CROSS-CORRELATION 

So far we have focused on the correlation function of the counts with themselves. Another useful cosmological probe 
is the correlation with the cosmic microwave background, primarily due to redshifting of photons as they propagate 
through the evolving potentials correlated to the source densities (the integrated Sachs- Wolfe effect; ISW) j43|. There 
are additional contributions to the correlation due to Doppler terms when CMB photons scatter during reioniza- 
tion p4| , lensing magnification |?9|, and also non-linear effects such as the Rees-Sciama effect in non- linearly evolving 
potentials, inhomogeneous reionization and the Sunyaev-Zel'dovich effect. We focus on the linear effects here, so that 
the fractional CMB temperature anisotropy sourced well after recombination is 

A T (n) » J V dr] eT T (rh ■ v + <tjj + </>) , (50) 

where r = t{vi) is the optical to scattering of CMB photons between t]a and 77, and — re~ T is the visibility. This ex- 
pression neglects small linear contributions from large-scale perturbations at recombination and quadrupole scattering 
that can easily be included in a full numerical linear calculation (as we do here). Cross-correlation with Eq. (|37]) gives 
a variety of terms, dominated by the cross-correlation of the source density and the CMB ISW, but in general with 
additional non-negligible correlations of the velocities, densities, magnification and potentials. The most important 
contributions have been calculated separately before, but doing a consistent linear analysis ensures that no relevant 
effects arc missed. Additional velocity and magnification terms are relatively most important when the dominant 
signal is small, i.e. for sources at high redshift (when the universe is close to matter-dominated) so that the direct 
ISW contribution from that redshift is small. 

Contributions to the cross-correlation are shown in Fig. ^. At high redshift, where the potentials are nearly constant, 
the intrinsic correlation is small and the signal is dominated by lensing unless 2 — 5s is small in which case magnification 
bias partly cancels the change in angular density of sources from convergence. At high redshift the electron density is 



also higher, giving a significant contribution from the CMB Doppler source [g9, 44 . At lower redshifts the correlation 
of the density and ISW dominates as expected. In addition to the late-time Doppler and ISW sources, there is also 
a correlation with the Sachs- Wolfe signal from recombination due to very large Hubble-scale modes; these correlate 
both with the density and velocity count sources on very large scales at the < 10% level (falling rapidly from 1 = 2). 

On all scales the redshift-distortion and other velocity contributions to the count part of the correlation are small, 
contrary to the conclusion of Ref. Q. Redshift distortions would increase the apparent density of sources over 
the peak of a matter overdensity, but lead to redshift-space underdensities in the tails even though the potentials 
contributing to the ISW are still large and the same sign there. Since the redshift-space density averages to zero, the 
overall correlation is close to zero. Mathematically, for a given / and sharp window at \ — X*i the oscillating j'/ (kx*) 
redshift-distortion term in Eq. ( [io|) is integrated over lnfc against a smooth source for the ISW [oc f dr)(<j) + ip)ji(kx)], 
giving close to zero; see Fig. [?]. However for numerical work the small few-percent correction on very large scales can 
easily be included. 

In addition to the correlation between the CMB temperature and counts, there is also some large-scale correlation 
with the E'-mode polarization generated by scattering at reionization. The quadrupole seen by an electron at reion- 
ization is generated by the Sachs- Wolfe effect on the electron's last scattering surface; since the Sachs- Wolfe effect 
is proportional to the gravitational potential and the horizon size at reionization is large (and hence the correlation 
length corresponds to I ~ 6), the polarization generated by Thomson scattering is correlated with the large-scale 
potentials at significantly lower redshift. For further discussion in the context of lensing see Ref. ^(J. For counts at 
z < 3 the correlation is in principle marginally detectable on the full sky at ~ 2<r, but falls rapidly on small scales 
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FIG. 6: Cross-correlation of the CMB temperature and number counts at z — 0.6, a z = 0.05,6 = 1, s = (left) and z = 
3, <7 Z = 0.2, b — 2, s = 0.42 (right) assuming no source evolution. Contributions are density-ISW (dashed black), magnification- 
ISW (dot-dashed red; absolute value), density-Doppler (solid blue), radial velocity gradient-ISW (dashed magenta; redshift 
distortions), time delay-ISW (thick cyan), total counts-CMB Sachs-Wolfe (thick magenta marked LSS; the CMB contribution 
is from recombination), and total (thick solid black). The left figure shows the absolute value of the contributions on a 
logarithmic scale, and additionally shows the contribution from terms involving non-integrated gravitational potentials (dotted 
green) and velocity-ISW (lower dot-dashed blue). The lensing contribution with s = is negative; the contribution at z — 3 is 
low here (right) because 5s — 2 = 0.1. 




fcMpc 

FIG. 7: Transfer functions at I — 20 for the ISW (thick solid black) and Doppler (dot-dashed black) contributions to the CMB 
temperature, and density (thin dashed red) and redshift-distortion (thin solid magenta) contributions to the counts source for a 
window function peaked at z = 0.6 with a z = 0.05. The corresponding contributions to the Ci are the integrals of the products 
over lnfc against the primordial power spectrum which is nearly flat. The Doppler and redshift-distortion contributions mostly 
cancel. CMB transfer functions are scaled by 5 x 10 . 
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FIG. 8: The correlation coefficient r = Cf A / ^Cf E Cf lA between counts (A) and CMB _E-mode polarization for various 
redshifts, assuming Gaussian redshift bins of width a z = 0.5, unit bias, and reionization redshift z IC = 11. 



and low redshift as shown in Fig. g. However the correlation is significantly larger than calculated in Ref. [[17) where 
only the much smaller low-redshift signal from re-scattering of the ISW signal was included. Combining counts at 
z < 3 with CMB lensing to reconstruct the higher redshift potentials, the correlation with polarization is in principle 
detectable at over 6a, and hence should be accounted for in any self-consistent full joint analysis. 

VII. CONCLUSIONS 

We have shown exactly how observable angular source densities as a function of redshift relate to the underlying 
physical densities, velocities and potentials. In linear theory we recovered various well-known terms, as well as new 
linear terms that are negligible for modes well inside the Hubble radius. The gauge-invariant observable results include 
in a consistent way the different terms that can be described in terms of gravitational lensing, source evolution, radial 
displacement and other velocity and general- relativistic effects. Our numerical code is available 6 , allowing rapid 
calculation of the auto and cross-spectra for multiple radial window functions, with cosmic shear, and with the CMB 
and the CMB lensing potential. 

We have not addressed in detail the more difficult question of how the number of sources relates to the underlying 
cosmological perturbations and background cosmology. Nor have we included non-linear effects. Several of these are 
in fact much more important than the small extra velocity and post-Newtonian terms that we derive here. Apart 
from non-linear evolution effects on the distribution of the source numbers (and velocities) themselves, there is also 
lensing of the perturbations |}0| and non-linearities in the redshift to real-space mapping [30| , both of which 
can have a significant effect on the power spectrum and correlation functions on small scales. 
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Appendix A: Relating number counts and line radiation 

Another important observable is the observed brightness from diffuse line emission (or absorption), especially 21-cm 
radiation from the spin- flip transition in hydrogen p2| . A full linear-theory treatment is complicated due to line-of- 
sight effects and is discussed in detail Rcf. ^(J. Here we make the connection with the differential number counts 
discussed in the main text. 

In their rest frame, and neglecting the finite line width, the sources with number density n s emit line radiation 
isotropically at frequency Eq. In proper time dr the number of photons added per volume within frequency dE is 

an s 5(E - E )dEdr, (Al) 

which defines a. We shall assume a is independent of spacetime position, and neglect scattering and absorption after 
emission. For a more detailed discussion see Ref. Q . 

We consider collecting the line photons by an observer with 4- velocity u\, equipped with a detector sensitive to 
photons in an energy range dE about E and with area dA that admits photons in solid angle d£l. The number of 
photons collected in proper time dt for this observer is 

dAf = f{E)E 2 dEdndAdt, (A2) 

in terms of the photon distribution function /. Emitting material around affine parameter A on the line of sight is at 
redshift z. It emits in its rest-frame at energy Eq and this will contribute to cLV" for E(l + z) = Eq. The (invariant) 
area of the wavefront associated with rays in the bundle dfi is dA = detl?(A)df2 at A, and the collecting area dA 
subtends a solid angle dCt in the rest-frame of the emitter at z. As the ray advances by dA, the wavefront sweeps out 
a volume dAit°fc Q dA and it adds to the photons that will be detected in dE and dt by an amount 

an s {l + z)dE5[E(l + z) - E ]k a u a s d\dA^-— — . (A3) 

47T 1 + Z 



It follows that 



cW = adtdE J dA ( k a J a S[E(l + z) - E^dA^- ) ■ (A4) 



The 5- function allows us to evaluate the integral simply. If we then compare with Eq. (A2) and use the reciprocity 
relation, dAdVL = dAd£!/(l + zf, we find 



Ef(E, n) =ak a J a 



dA 
dz~ 



n{n,z) 

a l^D> (A5) 



where E(l + z) = Eq and a = q/(47t_Eq) is another constant. It is straightforward to show that if we use Eqs (E 



and ( |23| ) in the right-hand side of this expression, we recover the result in Ref. |20 when line-of-sight scattering and 
self- absorption is neglected. 
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